Nonlocal-local multimode bifurcation in turbulence 
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It is shown that a mechanism of energy redistribution and dissipation by the inertial waves 
can be effectively utilized in isotropic turbulence at small Reynolds numbers. This mechanism 
totally suppresses the local interactions (cascades) in isotropic turbulence at the Taylor-scale based 
Reynolds number 7?^ < 10'^. This value of R\ (that can be considered as a bifurcation value at which 
the local regime emerges from the nonlocal one in isotropic turbulence) is in agreement with recent 
direct numerical simulations data. Applicability of this approach to channel flows is also briefly 
discussed. A theory of multimode bifurcations has been developed in order to explain anomalous 
(in comparison with the Landau-Hopf bifurcations) properties of the nonlocal-local bifurcation in 
isotropic turbulence. 

PACS numbers: 47.27.-i, 47.27.Gs 



INTRODUCTION 

Thanks to increase of resolution of direct numerical 
simulations (DNS) the significant problem of emergence 
of the local interactions in isotropic turbulence became 
available for a quantitative investigation (see, for in- 
stance, Refs. [H-Q and references therein). Using the 
DNS data it is shown in a recent paper Q that nonlocal 
interactions are completely dominating in low Reynolds 
number isotropic turbulence, as well as in a near dissi- 
pation range at moderate and large Reynolds numbers. 
The bottleneck phenomenon produced by these 

interactions, has an universal character. In correspond- 
ing near-dissipation range of scales 5]: rj < r < rjs the 
nonlocal interactions are completely dominating ones in 
isotropic turbulence (where r^s — 27ry and 77 is the dissi- 
pation or Kolmogorov scale). Another scale: Tc, is a 'sta- 
bility exchange'" scale [H. Namely, for scales r > Tc, the 
local (cascade) regime is stable and the nonlocal regime 
is unstable, whereas for scales r < Tc the local regime is 
unstable and the nonlocal regime is stable. The normal- 
ized scale rjs /rj is constant [sj whereas rdrj increases with 
Rx [1] ■ Therefore, there is a value of R\ at which Tc ~ rjs 
(see Figs. 1 and 2). This value of R\ can be considered 
as a bifurcation value of Reynolds number at which the 
local regime (still unstable) emerges from the nonlocal 
one. In present paper a phenomenological approach is 
suggested in order to compute bifurcation value of R\ . 

There is an overlapping of the local and nonlocal 
regimes in the range of scales Tc > r > r^g Q , Q . More- 
over, in the turbulent environment even the bifurcation 
scale Tjs might fluctuate around its mean-field value along 
with the dissipation scale ry (see 0] , , [11 ) • Hence, at 
the nonlocal-local bifurcation the supercritical and sub- 
critical secondary oscillations should co-exist. However, 
at the Landau-Hopf bifurcation this cannot be possible 
(see, for instance, Refs. HI)- In order to explain 

this 'anomalous' behavior of the nonlocal-local bifurca- 
tion we will develop a theory of multimode bifurcations. 
In this theory several modes (which frequencies are in cer- 



tain specific relationships) are exited simultaneously. At 
a bifurcation in a turbulent environment at least second 
harmonic should be taken into account. From dynamical 
point of view this situation is more complex than the one- 
mode (Landau-Hopf) bifurcation and corresponding dy- 
namical equations for the exited modes are multidimen- 
sional (in the one-mode Landau-Hopf bifurcation the dy- 
namical equation for the exited mode is two-dimensional 
in real variables). It is well known that behavior of dy- 
namical systems with dimension more than two is sig- 
nificantly different from the two-dimensional dynamical 
systems. In particular, the former equations can generate 
stochastic attractors. This gives to the multimodal bifur- 
cations considerable capacity in the phase space, which 
(unlike to the one-mode bifurcations) provide the mul- 
timode bifurcations a possibility to survive in turbulent 
environment. 



PHENOMENOLOGY 

A large vortex produces a Coriolis force field in the 
domain of space that it occupies, which is due to the 
rotational motion of the vortex. The comparatively slow 
motion of a large vortex affords time for the process of ra- 
diation of inertial waves to be realized by small-scale fluc- 
tuations in the Coriolis force field of this vortex [l^] , [13] ■ 
Such inertial waves transfer the energy being emitted 
by the fluctuations to the viscous layers on solid walls 
or to the viscous layers generated by the inertial waves 
on the effective boundaries of the large scale vortices 
14j . It is known that the energy brought by the inertial 
waves to the viscous layers dissipates effectively therein 
14| . The dissipation layer of inertial waves has 



mm 

a thickness of order 



/ (X (1^/17)1/2, 



(1) 



where v is the molecular viscosity, and f2 is a char- 
acteristic angular velocity of the large-scale vortex. To 
estimate the value of il, one can use £7 oc iT^u^jL. where 



2 



Mo is a typical large-scale velocity, and L is a typical scale 
of the large-scale vortices. Then 



I cx (y i/27ruo)^/^ 



(2) 



If we define large-scale Reynolds number = Luq/v, 
then we obtain 



/ cx 



2ttR 



1/2- 



(3) 



Or, using the Taylor-scale based Reynolds number Rx 



I 



2nRx 



(4) 



The dissipative layer of the inertial waves is the 'friction' 
layer of the corresponding large-scale vortices. 

Part of the energy entrained by the inertial waves dis- 
sipates in the viscous layers while another part, being 
reflected, returns to the volume occupied by the bulk of 
the vortex. Decay of turbulent fluctuation kinetic energy 
a result of its removal by inertial waves to the 'friction' 
layers can be described by equation 



dt 



1/2 



(5) 



where a is a dimensionless constant. If and L are 
considered constants (or slowly changing), then it follows 
from (5): 



(X exp - 



At 



(6) 



where A = a(j/il)^/^/L . This decay can be taken into 
account in the equation for the spectral tensor in the 
"external friction" approximation [12| : 



dt 



T^j{k,t)-\E, 



(7) 



where Tij characterizes the spectral energy transfer due 
to the nonlinear effects. The "external friction" reflects 
the interaction between inertial effects and the viscous 
('friction') layers, and it is a phenomenological substitute 
to the viscous term. Since the inertial term in the Navier- 
Stokes equation is quadratic in the velocity, then for small 
Eij the functional Tij is a homogeneous functional of 
Eij of order 3/2. Let us make the substitution 



E*j{k,t) ^ E,j{k,t)e 
Substituting (8) into (7) we obtain 



dE, 



dt 



Now, we make a replacement of the time 

= ^(l-e-^*/2). 
A 



(8) 



(9) 



(10) 
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FIG. 1: Scales rc/rj (solid line) [§] and rjs/ri 0] (dashed line) 
vs. R\. The arrow indicates the transitional value of the 
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Then it follows from (9) that 



dE^ 
dt' 



Tij (fc, t ). 



(11) 



Replacements (8) and (10) reduce the problem with "ex- 
ternal friction" to a problem without friction (the initial 
conditions are evidently identical). As t ^ oo , it follows 
from (10) that t 2/X, so that the whole evolution 
with "external friction", Eq. (7), is stacked in the in- 
terval (0, 2/A) of the evolution described by ideal Eq. 
(11). 

On the other hand, it is known (see, for instance, Refs. 
[l3|i[3l'0) that there is a certain characteristic time of 
an inviscid cascade process development 



tc 



Uo 



(12). 



Taking the preceding into account, an estimate can be 
made of the condition under which the cascade process 
cannot be developed successfully (because of the action 
of the " external friction" ) . This condition has the form 



2/A <t,~2—, 

Uo 



2L/a{uny/'^ < 2—. 

Uq 



(13) 



(14) 



Using estimate a ~ 10 [12l | we can rewrite condition 
(14) in the form 



Rx < 10^ 



(15) 



Thus, in isotropic turbulence with Rx < 10^ the cas- 
cade process (local interactions) cannot be developed suc- 
cessfully (cf. Ref. [20|). Figure 1 shows dependence of 
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R,=54 




FIG. 2: Three-dimensional spectra from a DNS performed in 
Pi for different Reynolds numbers. The solid curves in this 
figure correspond to the best fit by nonlocal interactions [^. 



the scales Tc/rj 8] and rjs/rj 5] on Rx. Intersection of 
the lines: rc/r/ vs. R\ and 7]s/'n vs. indicates the 
transitional value of R\ ~ 10^ (see Eq. (15) and Intro- 
duction) . Figure 2 illustrates the transition as it appears 
in the energy spectra of isotropic turbulence. 



SEPARATION 

One can also try to apply the above consideration 
to channel flows. The large-scale vortices acquire a 
metastable state in channel turbulence because of the 
action of inertial forces. These forces are equilibrated by 
viscous friction on the solid walls of the motion domain. 
The dissipative layer of the inertial waves on the solid 
walls is the 'friction' layer of the corresponding large- 
scale vortices. Therefore it is plausible to assume that 
I Up, where Wj, is the characteristic size of the 'criti- 
cal' layer [23|,[22[. With this assumption we obtain from 
Eq. (3) 



Vp oc R 



1/2 



(16) 



where y+ = j/p u,,/v, is so-called friction velocity, 
and i?, — L u^jv 2l| (here L is a characteristic size 
of the channel, and ~ uq). Eg. (16) is in agreement 
with the data presented in Ref. j21l]. This can be consid- 
ered as a positive indication of applicability of the above 
developed approach to the channel flows as well. How- 
ever, in order to estimate transitional (from nonlocal to 



local turbulence) value of Reynolds number for the chan- 
nel flows one needs in an analogy of the estimate (12) for 
such flows ufh. 



MULTIMODE BIFURCATIONS 

Existence of the overlapping range of scales: Tc > 
r > r]s, implies co-existence of the supercritical and sub- 
critical secondary oscillations at the nonlocal-local bi- 
furcation in turbulent media. At the standard Landau- 
Hopf bifurcation [l3,[lH such co-existence is impossible. 
Therefore, one needs in a generalization of the bifurca- 



tion theory (see Ref. [23|,[2j| for other examples of bifur- 
cations in turbulent and noisy media). The bifurcation 
problem for stationary solution of an evolutionary nonlin- 
ear equation in a Hilbert space generally can be reduced 
to investigation of instability of trivial solution (u = 0) 
of an equation 



du f 



N 

E 

n=2 



(17) 



where u is an element of some Hilbert space, L is a linear 
operator in this space, Fn{u) are nonlinear operators in 
this space such that _F„(Au) = A"-F„('u), and A'' is called 
order of nonlinearity. Generally the linear operator L is 
a nonself adjoint operator. Complete system of vectors 
in the Hilbert space contains eigenstates and associated 
eigenstates of this operator [2^ . Let us expand a solution 
of the equation (17) in series of this full system of the 
vectors 



i{t) = ^ ak{t)'ipk 



(18) 



fe=i 



Since physical fields can take only real values the spec- 
trum of operator L always contains pairs of complex- 
conjugate eigenvalues 25|: Lip = and Lip = jip, 



where 7 is complex conjugate to 7. 

Let us start from the one-mode case. Let to a pair 
of such eigenvalues: = 6 + zwq, and 71 (correspond- 
ing to the pair of eigenstates oi and 02) be located in a 
straight line which placed in a small vicinity of the imag- 
inary axis (i.e. 1 \S\). Rest of the eigenvalues of the 
operator L are located on the left side of this straight 
line and of the imaginary axis. The parameter 6 is con- 
trol parameter of the system {S ~ {Rl — Rl,c), and Rl,c 
is a bifurcation value of the large-scale Reynolds num- 
ber) . For small enough amplitude of the oscillations we 
can restrict ourselves by consideration of = 2. Let 71 
be simple eigenvalues [25[- Then, substituting the series 
(18) into (17) and multiplying on the eigenvectors of the 
operator L"*" (adjoint to L) in the Hilbert space we obtain 
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for ai equation: 

da I 
It 



71 a, 



k,m—l 



(19) 



bkm = i(f>l,F2{lpk,1pm))/i(t>l,1pl); 

where are eigenvectors of the operator L"*". We have 
taken into account that the vector 0i is orthogonal to aU 
tpk excluding ^pl . Equation for a2 —oT can be obtained 
using the complex conjugate operation. 

Let us introduce the variable Q — tot, where lu is fre- 
quency of the auto-oscillations. Let us expend a„(0) and 
LU in analytic series on S 



LU = ajQ + Suii + S''uj2 + 



(20) 



(21) 



Substituting these expansions into equations for a„ we 
obtain in the first order 



a[^\e) = ^iexp(ie) 



(22) 



Amplitude Ai can be obtained from next orders of the 
expansion. It can be shown that for n > 2 



e 



aW(e) = (co + cie-f... + Cfee'=)exp(7„ — ) (23) 

(the polynomial multiplier is related to the associated 
eigenvectors [2^). Because 3fJ{7„} < for n > 2, the 
an\Q) decay with time for n > 2. Further the terms 
decaying with time will not be taken into account. 
In the second order we obtain 



J (2) 

dal . (2) 
— la] ' = — 

de ^ ^„ 

n.m— 1 

(24) 

where b^^^^ = bkm{S = 0), /? = 1 for supercritical and 
(3 — —1 for subcritical cases. Analogous equation can 

(2) 

be obtained for the complex conjugate function by 
complex conjugation operation. 

If one substitutes a^^^ from (22) into right-hand side 
of Eq. (24), then terms ~ exp(iO) and ^ exp(— zO) will 
lead to resonance, which cannot be eliminated due to the 
necessary for this elimination condition 



(25) 



cannot be satisfied. Therefore the analytical series ex- 
pansions (20), (21) are not applicable in the one-mode 
case. 

Instead of the analytic expansions (20), (21) one can 
use non-analytical (Hopf) series on S^^^ 



(26) 



LO = UjQ + + 5lu2 + 



(27) 



In the first order we again obtain for a^^^ representation 
(22). But in the next order we now obtain equation 



da{ 



(2) 



-,(2) 



(1) 



^ E (28) 



uji da{ 

luq dQ ujQ 



n,m — l 



The resonance term can be now eliminated by choosing 
uji = 0. Analogous situation takes place for complex 

(2) 

conjugated equation for 03 . Then in the next order 



da 



(3) 



(3) _ ^2 da\ 



(1) 



de 



iP\ ' 1 

ojQ d@ loq 



2 

E 

n.m— 1 



(29) 



Since the last term in equation (28) gives terms ~ 
exp(z20) and ~ exp— (j28) into its solution the reso- 
nance terms can now be eliminated from the equation 
(29) and conditions of this elimination provide us Ai 
and LU2. This is _ff op/ bifurcation which gives dependence 
ai ^ ,51/2. 

If the original stationary solution has groups of sym- 
metry, then more than one pair of the eigenvalues can be 
located in a straight line parallel to the imaginary axis. 
Let to two pairs of such eigenvalues: — S + iujQ, 
and 72 = ^ -l- 2ict;o, 72, be located in a straight line which 
placed in a small vicinity of of the imaginary axis. Rest 
of the eigenvalues of the operator L are located on the 
left side of this straight line and of the imaginary axis. It 
will be clear from further considerations that fulfillment 
of these conditions with accuracy 0{S^) is enough for ex- 
istence of the auto-oscillations in the system. Let 71 and 
72 be simple eigenvalues ^5]. Then, substituting the se- 
ries (18) into (17) and multiplying on the eigenvectors of 
the operator L+ (adjoint to L) in the Hilbert space we 
obtain for ai and 02 following equations: 



daj 
~dt 



(30) 



k.m—l 



bjkrn = {(/)], F2{lpk, 1pm)) /i(f>j,i^]); 

where (pj are eigenvectors of the operator L+. We have 
taken into account that the vector (pj is orthogonal to all 
ipk excluding ipi, and 4>2 is orthogonal to all ipk exclud- 
ing ip2- Equations for as = oT and for a4 = 02 can be 
obtained using the complex conjugate operation. 

Let us now expend a„(0) and lo in the analytic series 
(20), (21). Substituting these expansions into equations 
for a„ we obtain in the first order 



z$.^)(e) = A,exp(z5Me) 



(31) 
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Amplitude Aj can be obtained from next orders of the 
expansion. In the second order we obtain 

dQ loq ■' ujQ dQ uq ^ 



i- E ^y^'^ (32) 

n,m— 1 

If one substitutes o'p taken from (31) into right hand side 
of equations (32), then terms exp(i0) and ~ exp(i20) 
will lead to resonances for j = 1 and for j = 2 respec- 
tively. If we write conditions of suppression of these res- 
onances 

(z2wi - /3)wo 
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FIG. 3: Dependence of v! on i?J for a grid-turbulence exper- 
iment l26|. The straight line indicates Eq. (36). 



fc = fci + ^k, = -5R{(6l?3 + &f32)4i\} + mb^i2\ + bfl\} 



we obtain 



l^ll =-1— Wl, 



3fci^^| /9fc2 
4fc2 f3 



16fc| 



(33) 



(34) 



We have found only modulus of Ai , but we can let Ai be 
real due to initial phase of the two- mode auto-oscillations 
can be arbitrary chosen. 

Thus one can see that for the multimode bifurcation 
the amplitude of the secondary oscillations is propor- 
tional to the control parameter, S, while for one-mode 
(Landau-Hopf ) bifurcation this amplitude is proportional 
to Vs. It is clear that the anomalous dependence will 
also take place in, for example, more complex situation 
of three complex conjugated pairs with frequencies luq, 
(n — l)uJo, and nujQ (where n is an arbitrary integer). 
One can construct other combinations of the frequencies 
leading to the anomalous dependence. 

In our case 5 ~ {Rl—Rl,c) (where Rl,c is a bifurcation 
value of the large-scale Reynolds number). Let us denote 
for longitudinal fluctuations of the velocity field v a rms- 
value: u' = (u^)^/^. Then for the multimode bifurcation 



u — + C{Rl — Rl,c 



u'^u, + ^C{Rl - i?2 J -t- 



(35) 



(36) 



(where C is certain constant). Hence, for the multi- 
mode bifurcation u' is a linear function of i?^ in a vicin- 
ity of the bifurcation value i?A,c (fo r the standard one- 
mode Landau-Hopf bifurcation [l3|,[lll we have u' = 
+ C^Rl — Rl.c: instead of the Eq. (35)). Figure 
3 shows dependence v! on for the data obtained in a 
grid-turbulence laboratory experiment [3] ■ The straight 
line is drawn in this figure in order to indicate agreement 
with the Eq. (36) (i.e. with the multimode bifurcation 
prediction). 



STABILITY 

It follows from Eq. (33) that its right hand side should 
be positive. It can be obtained for any sign of /? if one 
chooses an appropriate sign in (34). Therefore, the su- 
percritical and subcritical secondary oscillations co-exist 
at the multimode bifurcation (while at the standard one- 
mode Landau-Hopf bifurcations this cannot be possible). 
This co-existence makes the question about stability of 
the secondary regimes rather important. 

Let us consider the evolutionary equation linearized on 
the secondary regime 



(t^o 



(/3(S -f 7jt^o)aj 



fc,m— 1 



(37) 



where (0) are the first order (on the parameter S) 
modes of the considered secondary oscillations. The Flo- 
quet representation for evolutionary operator of these 
equations is [27i| 



1/(6) = Q(e)exp(Ge). 



(38) 
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Stability or instability of the secondary oscillation is de- 
termined by spectrum of the operator G. We use small 
parameter 5 to find when the spectrum of the operator 
G has its points located in the right semi-plane. Let us 
expend operator G 



I thank T. Nakano, D. Fukayama, T. Gotoh, and K. 
R. Sreenivasan for sharing their data and discussions. 



where 



G — Go ^" (5Gi 



Go — i diag{l, 2} 



(39) 



(40) 



Eigenvalues of the operator Gg are: 

■^2 



(0) 



i and 
(1,0), 



(72 ' = 2i. Corresponding eigenvectors are: xi 
X2 = (0, 1). To find operator Gi we use the Krein repre- 
sentation [271 



(A-^)e(^-^)" 

1 _ g(M-A)T 



x{Go~ \i)-^Ri{Go- fiiy^drdXdfi (41) 

where To = 2TrujQ^, 

Ri = ujQ^diag{f3 — iu)i, P — 2iuJi} 

and Tq is a contour enveloping the spectrum of the oper- 
ator Go . Harmonic perturbations of the right hand side 
of the evolutionary equation do not give a contribution 
to the 3?{tT} ~ 6. If we now expand spectrum of the 
operator G: a = cr'"-' + 6a^^^ + then we obtain 



rW = (y,Gix) 



(42) 



where x is eigenvalue of the operator Go, y is eigenvalue 
of the operator conjugate to Go and (,) means scalar 
multiplication. Then, substituting Gi from (41) into (42) 
and calculating the integrals we obtain 



5R{(7} 



(43) 



Therefore, for small 5 the subcritical secondary oscil- 
lations (/3 = —1) are unstable (3?{ct} > 0), whereas 
supercritical {(3 — 1) secondary oscillations are stable 
{lk{a} < 0). 

While one-mode equation (19) is two-dimensional (for 
real variables) the multimode equations are multidimen- 
sional: two-mode equation has dimensionality 4, and the 
three-mode equation for the modes with ujq, (n — 1)luo 
and uluq has dimensionality 6. It is well known that non- 
linear equations with dimensionality larger than 2 and 
equations with dimensionality smaller (or equal) than 2 
have very different properties. In particular, the former 
equations can generate stochastic attractors. These at- 
tractors have considerable capacity in the phase space 
that allows them to survive in the turbulent environment. 
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